LDOS modulations in cuprate superconductors with competing AF order: the 

temperature effect 
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Based upon a phenomenological t — t' — U — V model and using Bogoliubov-de Gennes equations, 
we found that near the optimal doping S — 0.15 at low temperature (T), only the pure d-wave 
superconductivity (dSC) prevails and the antiferromagnetic (AF) order is completely suppressed. 
However, at higher T the AF order with stripe modulation and the accompanying charge density 
wave (CDW) emerge, and they could exist even above the superconducting transition temperature. 
This implies that the existence of the CDW depends critically on the presence of the AF order, not 
so much on the dSC. The LDOS (local density of states) image indicates that the stripe modulation 
has an energy independent spacing of 5a /4a spreading over a 24a x 48a lattice, corresponding to 
an average periodicity 4.8a. This result may be relevant to the recent STM experiment [Vershinin 
etal., Science 303, 1995 (2004)]. 

PACS numbers: 74.25.Jb, 74.20.-z, 74.50.+r 



Recently, STM measurement by Vershinin et al. |l| on 
slightly underdopcd BSCCO indicated that the electronic 
states at low energies and at a temperature T higher than 
the superconducting transition temperature (T c ) in the 
pscudogap region exhibit an energy-independent spatial 
modulation which resulting to a checkerboard pattern 
with incommensurate periodicity 4.7a ±0.2 (a is the lat- 
tice constant). At very low temperature, however, no 
such pattern has been detected in agreement with 
previous measurements 0, 0. So far the reasons why 
this STM pattern at low temperature looks so different 
from that at high temperature, and why the noninteger 
incommensurate periodicity could occurred in a under- 
line Cu02 lattice are outstanding questions which have 
not been addressed in the existing literatures. In the 
present paper we are trying to understand these issues 
by adopting the idea of the d-wave superconductivity 
(dSC) with a competing antiferromagnetic (AF) order for 
cuprate superconductors 0, 0, S S > an d to exam- 
ine the formation of the AF order and the accompanying 
charge density wave (CDW) at finite temperature. Em- 
ploying the phenomenological t — t' — U — V model with 
proper chosen parameters, we show that at low temper- 
ature only dSC prevails in our system and the AF order 
is completely suppressed. The LDOS (local density of 
states) image is featureless. At higher temperature, it is 
found that the AF order with stripe modulation, which 
is referred also as the spin density wave (SDW), and the 
accompanying charge density wave (CDW) may show up 
and they could even persist at temperatures above the 
BCS superconducting transition temperature T^ cs . In 
the presence of SDW, we show the LDOS image to have 
an energy-independent stripe modulation with spacing 
of 5a /4a spreading over a 24a x 48a lattice. Accord- 
ing to the Fourier analysis of the LDOS image, an av- 
erage periodicity 4.8a could be assigned for the stripe 
modulation. If the doubly degenerate states of a;-and y- 



oriented stripes have the same probability to appear in 
the time interval when the experiments is performed, the 
combined LDOS image should have a checkerboard pat- 
tern of 4.8a x 4.8a structure. This result is in agreement 
with the experiment by Vershinin et al. Q. In addition, 
the temperature dependent behavior of the normalized 
LDOS, exhibits the "pseudogap"-like characteristics 0| 
originated in the emergence of the SDW at higher T. 

To model these observed phenomena, we employ an 
effective mean-field t — t' — U — V Hamiltonian by as- 
suming that the on-site repulsion U is responsible for the 
competing antiferromagnetism and the nearest-neighbor 
attraction V causes the d-wave superconducting pairing 
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where iy is the hopping integral, /i is the chemical po- 
tential, and Ay = -|-( Ci T c ji — c U c jt) ^ s the spin-singlet 
d-wave bond order parameter. The Hamiltonian above 
shall be diagonalized by using Bogoliubov-de Gennes' 
(BdG) equations, 
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where Hij a = — + (U(ni a ) — Here, we choose 

the nearest- neighbor hopping (ty) = t = 1 and the next- 
nearest-neighbor hopping (iy) = t' = —0.25 to match 
the curvature of the Fermi surface for most cuprate su- 
perconductors [l2^ . The exact diagonalization method to 
sclf-consistcntly solve BdG equations with the periodic 
boundary conditions is employed to get the N positive 
eigenvalues (E n ) with eigenvectors (u-i, vfA and N nega- 



tive eigenvalues (E n ) with eigenvectors (- 



<). The 
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FIG. 1: The maximum value of the staggered magnetization 
(Mi)max as a function of U at zero temperature, where Mi = 

(-lAnit-mJ. 



self-consistent conditions are 
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1) is Fermi-Dirac dis- 



tribution function. Since the calculation is performed 
near the optimally doped regime, the filling factor, n/ = 
^2i<r{ c ia c ic)/NxNy, is fixed at 0.85, i.e., the hole doping 
S = 0.15. Each time when the on-site repulsion U or the 
temperature is varied, the chemical potential \i needs to 
be adjusted . 

In Fig. 1, we plot the maximum value of the staggered 
magnetization (Mi) max as a function of U with V = 1.0 
at T — 0. The value of U measure the strength of the 
competing AF order in the background of the dSC. In 
region I, where U > U c ~ 2.45, the AF order with stripe 
modulation or the SDW shows up and it has a period 8a 
oriented along the x- or y-axis. The corresponding modu- 
lations in dSC and accompanying CDW have a period 4a. 
All these have been discussed in one of our earliar works 
[l3j . In regions II and III, where U < U c , the AF order is 
completely suppressed by the dSC and (Mi) ma:E = 0, and 
the system is in the pure dSC state. The transition be- 
tween the stripe-modulated dSC/SDW/CDW coexisting 
state and the pure dSC state is discontinuous and of the 
first order. The difference between region II and region 
III is that in the presence of an applied magnetic field, 
dSC/SDW/CDW with stripe modulations are induced in 
region II while in region III two dimensional modulations 
are induced around the vortex cores. Detailed study of 
this issue has been given previously ^J] and will not be 
discussed here. In the following, we choose U = 2.44 
(region II) which is slightly smaller than the critical U c 
and U = 2.39 (region III) to examine their temperature 



FIG. 2: (a) and (b) show the temperature dependence of the 
dSC (open circle) and AF order (solid circle), respectively, 
for U = 2.44 and U = 2.39. The value of the dSC order 
parameter Ai = |(Ai_|_^ + Ai_^ — Ai + y — Ai_y) is measured 
in the unit of t and Mi — (M{) max . 



effects. 

Fig. 2 shows the temperature dependences of the 
dSC order (open circle) and the staggered magnetization 
(Mi) max (solid circle) for (a) U = 2.44 and (b) U = 2.39. 
From Fig. 2(a), it is straightforward to see that only 
dSC exists and the AF order is completely suppressed at 
low T. When T > 0.06t, the stripe-modulated AF order 
and the accompanying CDW start to emerge, and both 
of them could persist above the BCS transition tempera- 
ture = 0.14i. This result implies that the existence 
of the CDW depends critically on the presence of the 
AF order, not so much on the dSC. At T > T C BC5 , the 
staggered magnetization decreases rapidly to zero at the 
Neel's temperature TV = 0.155t. We also found that the 
stripe modulation associated with the SDW has mixed 
periods 10a/8a, and those associated with the dSC and 
the accompanying CDW have a mixed periods 5a/4a ori- 
ented either along the x-axis or the y-axis. Within the 
temperature range in which the AF order appears, these 
periods seem to be temperature insensitive. In Fig. 2(b), 
with U = 2.39, the AF order is completely suppressed 
and our system is in the state of pure dSC at all temper- 
atures. In both of these cases, the dSC order parameter 
as a function of T appears to have the BCS-like behav- 
ior. It is interesting to note that the T^ cs in the case 
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FIG. 3: Temperature dependence of the normalized LDOS 
(pi(E)t/pi(E) Tn , here Tat = 0.155t). The LDOS at T N as a 
function of E is represented by the dashed line. This repre- 
sentative set of spectra was shifted vertically for clarity. The 
wavevectors in first Brillouin Zone are M x x M y = 24 x 12. 



for U = 2.39 [Fig. 2(b)] is slightly larger than the one 
for U = 2.44 [Fig. 2(a)]. This is because the appear- 
ance of SDW in Fig. 2(a) at higher T also suppresses the 
dSC. The reason why the AF order is suppressed at low 
temperature and emerges at higher temperature could be 
understood from the phase diagram for U = 2.44 near the 
optimal doping which is shown in the inset of Fig. 2(b). 
The feature of this phase diagram is similar to what has 
been obtained for dSC with competing d-density wave 
(dDW) order 0. 
The LDOS is defined as 
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where pi(E)x is proportional to the local differential tun- 
neling conductance as measured by STM experiment, and 
the summation is averaged over M x x M y wavevectors in 
first Brillouin Zone. In order to investigate the contin- 
uous evolution of the quasiparticlc spectrum, the LDOS 



for U = 2.44 from low T to high T arc calculated. The 
normalized LDOS which is defined as p\(E)t / pi{E)r N , 
arc presented in Fig. 3 as a function of energy E 
and at temperatures ranging from T/Tn = 0.13 up to 
T/T N = 1.0. It should be noticed that T h c ulk (< T c scs ) 
labeled in Fig. 3 corresponds to the temperature where 
the coherent peaks of the dSC flatten out, and this sig- 
nature has been used to estimate the " superconductivity 
transition temperature" in many experiments. As it is 
indicated in Fig. 3 that the apparent "gap" is roughly 
a constant up to T h c ulk and progressively becomes slight 
larger up to TV, in consistent with the experimental mea- 
surements |xd . As far as the temperature behavior of the 
normalized LDOS is concerned, the essential feature of 
Fig. 3 displays the "pscudogap"-like characteristics. In 
the temperature range of T% ulk < T < Tjy, our system 
appears to be in the "pseudogap" region according to 
the characterization from STM experiments [Iollll| even 
though the effect due to the phase fluctuations on the 
dSC order parameter [^.ITsT] has not been taken into ac- 
count. In this region, the Fermi surface in our theory is 
everywhere gaped while the ARPES experiment [Tij in- 
dicates that the gaps occur only near (±7r, 0) and (0, ±7r). 
This difficulty so far has not been solved in the existing 
literatures. 

In order to compare with the STM experiment, we cal- 
culate the LDOS images between the energy 0.0 and 0.4 
with the increment 0.01. At the temperature T/Tn = 
0.94, which is larger than Tf cs and in the "pseudo- 
gap" region, we found that the LDOS images have stripe 
modulations with mixed periods of 5a/4a which is en- 
ergy independent. The periods of stripes spread into 
5a — 5a — 4a — 5a — 5a pattern oriented along either the x- 
axis or the y-axis on a 24 x 48 lattice. The system seems 
trying to establish a periodicity that is incommensurate 
with the underline lattice, but fails to do so because the 
calculation is performed in a discrete and finite lattice, 
not in a continuum. In Fig. 4(a), we show such an image 
at E = 0.1. The pattern in Fig. 4(a) still remains even 
as the temperature drops to T — 0.08t [Fig. 2(a)]. As 
it will be shown below that an average periodicity 4.8a 
for the stripe modulation could be assigned in this case. 
The x- and y-oriented stripe modulations are degenerate 
in energy and it is thus possible to use the minima of a 
double- well potential to represent the states of the x- and 
y- oriented stripes. Since the modulations discussed here 
are weak perturbations to the host system, we do not 
expect that the barrier in the mid of the double-well po- 
tential to be large and the tunneling between these two 
degenerate states could be achieved easily by quantum 
mechanics and thermal activation. When both x- and 
y-oriented stripe-modulations show up in the time du- 
ration when the experiment is performed, checkerboard 
patterns should be observed. 

In order to get a further understanding of stripe mod- 
ulation, let us investigate the Fourier transform of the 
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FIG. 4: (a) The LDOS image at E = 0.1 and T = 0.94Tjv. (b) 
The energy evolutions of the Fourier amplitude at T = 0.94Tjv 
(solid circle) and T = 0.86Tjv (open circle) of the LDOS image 
at q x = 0.417(2.), and the LDOS at T = 0.94Tjv. The insert 
in (b) is the Fourier amplitude of the LDOS image in (a). 
The wavevectors which have been used to calculate the LDOS 
image are M x x M y — 12 x 12 in first Brillouin Zone. 

LDOS image 

Pi(E)t = J— ^eip(iqTi) -pi(E) T ■ (5) 

The Fourier transform of the LDOS image in Fig. (4a) 
as a function q x is displayed in the insert of Fig. 4(b). It 
shows a very sharp peak around q x = 0.41 7( J) = jf^, 
corresponding to an average "periodicity" 4.8a in real 
space. At T = 0.94T W > Tf cs and in Fig. 4 (b), the 
curve represents the LDOS of a function of energy while 
the solid/open circle shows the energy evolution of the 
Fourier amplitude of the LDOS image at q x = 0.417(^), 
which arc normalized by their values at E = 0.4. Start- 
ing from the higher energy side, the solid circle first rises 
to a maximum at E = 0.27 and then dips as E ap- 
proaches to zero. This energy dependence is predicted 
for the case where we only have SDW/CDW without 
dSC in presence. When the temperature is lowered to 
T = 0.86TV < T® cs where dSC/SDW/CDW are all in 
coexistence, a result represented by the open circle is ob- 
tained for the Fourier amplitude of the LDOS image. If 
one introduces an additional parameter to consider the 
effect of phase fluctuations on the dSC order parameter, 
the dip in the low energy side of the open circle might be 
removed and flattened after reaching the maximum |20j 
and the obtained result could agree better with the STM 
measurement Q. A simplified procedure of doing this 
has been previously proposed |lSj and we are not going 
to repeat this again in the present paper. It also neces- 
sary to mention here that a different work has also been 
done by Chen et al. [2(j for rather underdoped sample ( 
S = 0.125) to address a similar issue. Their periodicity 
4a in the LDOS modulation shows up at all temperatures 
below the BCS transition temperature. 

In conclusion, we have investigated the interplay be- 
tween the dSC and the competing AF/CDW orders at 
finite temperature for samples close to the S = 0.15 dop- 
ing. With proper chosen parameters, we show that the 



AF order could be completely suppressed by the dSC at 
low T, and the LDOS image is featureless. At higher 
T, the AF and the accompanying CDW orders start to 
emerge and they persists even to temperatures higher 
than T C BCS . This finding implies that the presence of 
CDW depends on the existence of AF order, not so much 
on the dSC. When the AF order emerges, stripe like mod- 
ulations appear in all AF/dSC/CDW orders and also in 
the LDOS image. It is shown that the normalized LDOS 
from low T to high T exhibits the "pseudogap"-like be- 
havior. The Fourier transform of the LDOS image is 
peaked around q = (27r/4.8a,0) indicating a periodicity 
~ 4.8a for y-oriented stripes. Moreover the energy in- 
dependent checkerboard pattern which was observed by 
Vershinin et al. Q] at high T but not at low T could 
be understood in terms of the present theory if both the 
x- and j/-oriented stripes contribute in the time interval 
when the experiment is performed. Finally we point out 
that for a different sample which may have a U slightly 
larger than 2.44 (see Fig. 1 in region I) or a hole-doping 
level slightly less than 0.15, stripe modulations with pe- 
riod 4a would show up at T = 0. This result should 
be consistent with the energy-independent checkerboard 
patterns observed in other STM experiments [2~D |. 
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